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ABSTRACT 

The formation of a circumplanetary disk and accretion of angular momentum onto a 
protoplanetary system are investigated using three-dimensional hydrodynamical sim- 
ulations. The local region around a protoplanet in a protoplanetary disk is considered 
with sufficient spatial resolution: the region from outside the Hill sphere to the Jovian 
radius is covered by the nested-grid method. To investigate the thermal effects of the 
circumplanetary disk, various equations of state are adopted. Large thermal energy 
around the protoplanet slightly changes the structure of the circumplanetary disk. 
Compared with a model adopting an isothermal equation of state, in a model with 
an adiabatic equation of state, the protoplanet's gas envelope extends farther, and a 
slightly thick disk appears near the protoplanet. However, different equations of state 
do not affect the acquisition process of angular momentum for the protoplanetary 
system. Thus, the specific angular momentum acquired by the system is fitted as a 
function only of the protoplanet's mass. A large fraction of the total angular momen- 
tum contributes to the formation of the circumplanetary disk. The disk forms only in 
a compact region in very close proximity to the protoplanet. Adapting the results to 
the solar system, the proto- Jupiter and Saturn have compact disks in the region of 
r < 21rjup (r < 0.028 m.jup) and r < 66rsat {r < 0.061 rH,Sat), respectively, where 
fjup (?'H,jup) and rsat (rn.Sat) are the Jovian and Saturnian (Hill) radius, respectively. 
The surface density has a peak in these regions due to the balance between centrifugal 
force and gravity of the protoplanet. The size of these disks corresponds well to the 
outermost orbit of regular satellites around Jupiter and Saturn. Regular satellites may 
form in such compact disks around proto-gas giant planets. 
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1 INTRODUCTION 



Since the first detection by iMavor fc Queloj l|l995h . over 
300 extrasolar planets (or exoplanets) have been observed. 
Almost all exoplanets are believed to be gas giant planets, 
like Jupiter and Saturn in the solar system, because massive 
planets are observed preferentially. Although these planets 
are believed to form in the disk surrounding the central star 
(i.e., the circumstellar disk or protoplanetary disk), their 
formation process is not yet fully understood. Under spheri- 
cal symmetry, including the radiative effect, t he evolution of 
gas gian t planets has been invest igated ( e.g 



Mizuno et al 



1978; M izunclligSOl: IStevensonl 198 2 : Bodc nheinier fc Pollack 
1989 : IPoUack et al lll996l : llkoma et al. Il2000l . I2OOI''). These 



studies showed rapid gas accretion onto the protoplanet af- 
ter the formation of a solid core with 5 — 20 . The gas 
giant planet acquires a large fraction of its mass in this gas 
accretion phase. Since gas flows into the gravitational sphere 



(i.e., Hill sphere) of the protoplanet with an angular momen- 
tum, a circumplanetary disk forms around the protoplanet 
in this phase. Then, satellites may form in this circumplan- 
etary disk. Regular satellites around a gas giant planet are 
considered to form according to a scenario similar to the 
formation of earth-like planets in the solar system: the pro- 
tosatellite forms by the accumulation through mutual col- 
lision of the satellitesimals that forms after the dust grains 
sink towar ds the equatorial plan e of the circumplanetary 
disk (e.g., Istevenson et al. Ill986l ). Since regular satellites 
are considered to be a by-product of gas-planet formation, 
they provide an important clue for understanding it. Thus, 
it is necessary to understand the formation of the circum- 
planetary disk, which is the site of the satellite formation in 
the framework of gas planet formation. However, the forma- 
tion and acquisition process of angular momentum in the 
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circumplanetary disk cannot be understood by spherically 
symmetric calculations. 

In two-dimensional simulations, protopla net evolution 
authors 



19871: iLubow et al. Ill99£ 


: lKlevlll999l: iD'Angelo et al. II2OO2I: 


Tanigawa & Watanabd 


2OO2I: Dobbs-Dixon et al. 120071") 



who showed that it is possible for a circumplanetary disk 
to form around a protoplanet. They also found prograde 
rotation of the protoplanetary system (protoplanet + cir- 
cumplanetary disk) against the orbital motion of the planet 
in the circumstellar disk. Recently, the evolution of a proto- 
planetary system was calculated in three dimensions (e.g. , 
Miyoshi ct al. 1999; Kley et al. 2001; Bate ct ah] l2003l : 
D'Angelo et ah 1120031 : iMachida et al. II2008I '). These studies 
showed that the gas flow pattern in three dimensions is qual- 
itatively different from that in two dimensions: gas flows 
into the protoplanet system only in the vertical direction in 
three-dimensional simulations. Thus, three-dimensional sim- 
ulations are necessary for investigating the gas flow around 
the protoplanet, angular momentum of the protoplanet sys- 
tem, and circumplanetary disk formation. Moreover, the cir- 
cumplanetary disk should be resolved with sufficient spatial 
resolution, because regular satellites are distributed only in 
very close proximity to the planet; the regular satellites of 
giant planets in the solar system are distributed only in the 
region of r ^ 60r-p, where Vp is the radius of the planet. On 
the other hand, the radius of the gravitational sphere (or Hill 
radius) of the planet extends up to ru ~ 1000 r-p. For exam- 
ple, the Hill radii for Jupiter and Saturn are th.j = 743 rp 
and rH,s = 1083 rp, respectively. Thus, to investigate a pro- 
toplanetary system, we must resolve not only the planetary 
radius ~ rp, but also a sufficiently remote region from the 
Hill sphere (r ^ rn, or r S> 1000 rp), requiring the reso- 
lution of spatial scales differing by a factor of more than 
1000. 

Recently, iD'Angelo et al. I (120031 ) and iMachida et al. I 

l|2008l ) used the nested-grid method to r esolve a wide spatial 
scale. Although lo'Angelo et al~l (|2003l ) calculated a proto- 
planetary system in three dimensions with sufficient spatial 
resolution, they could not investigate the region near the 
protoplane t because they adopted a sink cell in the region of 
r < O.lrn. IMachida et al. I (|2008l ) calculated the evolution 
of a protoplanetary system without a sink, and estimated 
the angular momentum accreted onto the system. However, 
they adopted the isothermal approximation. This approxi- 
mation is considered to be valid in the circumstellar disk 
and areas remote from the protoplanet because these re- 
gions are optically thin due to their lower gas density. How- 
ever, this approximation is not valid in the vicinity of the 
protoplanet, because the gas temperatu re increases adiabat- 
ically as the gas density b ecomes high (jMizuno et al.|[T978l : 
iLunine fc Stevensonlll982l ). 

In this study, the thermal evolutio n of a protoplane - 
tary system is modelled according to Mizuno et al.l (| 19781 ). 
and implemented into a three-dimensional nested-grid code. 
Using such a code, the system is resolved from a sufficient 
distance from the Hill sphere (r > 10 rn) to the Jovian ra- 
dius (r = 0.8 rp), and the angular momentum flowing into 
the protoplanetary system and circumplanetary disk are in- 
vestigated. In addition, the spatial resolution in this calcu - 
lati on is much higher than those in ID'Angelo et al. I (|2003l ) 
and IMachida et al. I (|2008l ) . The calculations show that the 



angular momentum acquired by the system is independent 
of the thermal evolution around the protoplanet after the 
protoplanet grows to M > O.lMj. A thin circumplanetary 
disk appears only in a compact region near the protoplanet 
(r < 50rp). The centrifugal radius derived from the angular 
momentum of the protoplanetary system corresponds well to 
the region in which the regular satellites of present gaseous 
planets are distributed. In addition, the surface density has 
a peak in this region. 

The structure of the paper is as follows. The frameworks 
of our models are given in §2, and the numerical method is 
described in §3. The numerical results are presented in §4. 
The discussion of circumplanetary disk and satellite forma- 
tion is presented in §5. 



2 MODEL 

2.1 Basic Equations 

In this study, a local region around a protoplanet 
is considered using the sh earing sheet model (e.g., 
iGoldreich fc Lvnden-BeUll 19651 ). in which the self-gravity of 
the disk is ignored. In addition, no physical viscosity is in- 
cluded, and the numerical viscosity can be ignored because 
it is sufficiently small. Thus, an inviscid gas disk model is 
adopted here. The orbit of the protoplanet is assumed to 
be circular on the equatorial plane of the circumstellar disk. 
Local rotating Cartesian coordinates with the origin at the 
protoplanet are set up, in which the x-, y-, and z-axis are the 
radial, azimuthal, and vertical direction of the disk, respec- 
tively. The equations of hydrodynamics without self-gravity 
are solved: 



|+V.(p«) = 0, 

dv , „, 
- + (..V). 



— VP- V<E>eff 

p 



(1) 



(2) 



where p, v, P, $off , and i7p are the gas density, velocity, gas 
pressure, effective potential, and Keplerian angular veloc- 
ity of the protoplanet, respectively. For the gas pressure, a 
barotropic equation of state is adopted (for details see ii2.2|l . 
In the above equations, the curvature terms are neglected. 
The Keplerian angular velocity of the protoplanet is given 
by 



f^p 



GMc 



1/2 



(3) 



where G, Mc, and Op are the gravitational constant, mass 
of the central star, and orbital radius of the protoplanet, 
respectively. The effective potential "l>cff is given by 



ii(3x^ 



GMp 



(4) 



where Mp and r are the mass of the protopla net, and the dis- 
tance from the centre of the protoplanet (e.g.. lMivoshi et al. I 
I1999I ). The first term is composed of the gravitational po- 
tential of the central star and the centrifugal potential, and 
higher orders m x, y and z are neglected. The second term 
is the gravitational potential of the protoplanet. Using the 
Hill radius 

1/3 



rn = 



( Mp 



(5) 
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equation Q can be rewritten as 



fin 



r 



(6) 



2.2 Equation of State 



At a fixed orbital radius of tlie protoplanet, tlie circumstel- 
lar di sk lias an almost constant temperature l|Havaslii et al. I 
1 19851 ). while the gas around the protoplanet (i.e., the gas en- 



(iMizuno et al.lll97S 


: lMizundll980l: iBodenheimer & Pollack 


1 19861: Pollack et al. 


Il996l: Ikoma et al. II2OO0I). Mizuno et al. 



around the protoplanet, on the assumption that the enve- 
lope is spherically symmetric, and in hydrostatic equilib- 
rium. They also investigated the thermal evolution of the 
envelope, parameterising the dust grain opacity, and deter- 
mined the boundary between the isothermal and adiabatic 
regions. 

To correctly estimate thermal evolution around the pro- 
toplanet, we need to solve the radiation hydrodynamics. 
However, this has a hug e computational cost . Thus, in this 
study, using Figure 2 of iMizuno et al.l (|l978l ). thermal evo- 
lution around the protoplanet is modelled as a function of 
the gas density (i.e., barotropic equation of state) : 



2 

Cs,OP 



tanh 



P 



+ ftp^ tanh 



(7) 



where Cs is the sound speed, 7 is the adiabatic index (7 — 
1.4), and the adiabatic constant k is defined as 



2 1 — » 



(8) 



where pcri is the critical density: the gas behaves isother- 
mally in the region of p < pcri, and adiabatically in the re- 
gion of p > Pcri. In this study, pcri = 00 (isothermal model), 
lOpco, lOOpco, and lOOOpc.o (adiabatic models) are adopted, 
where pc,o is the initial density on the equatorial plane. The 
hyperbolic tangent (tanh) function is used to smoothly con- 
nect the first (isothermal) and second (adiabatic) terms in 
equation (O . The thermal evolution for different pcri is plot- 
ted against the gas density in Figure [T] in which the gas 
temperature is constant in the isothermal model (pcri = 00), 
while it increases gradually from the initial value at p ~ pcri 
in the adiabatic models (pcri = lO pc.o, lOOpc.o, and lOOO pc.o). 

In the standard disk model (|Havashi et al. Ill985h . the 
density and temperature at Jovian orbit are pc,o = 1.5 x 
10~^^ g cm~^ and To = 123 K, respectively. Thus, in a model 
with Pcri = lOpco (the highest temperature model), the gas 
behaves isothermally when p ^ 1.5 x 10"^" gcm~^, while it 
behaves adiabatically when p ^ 1.5 x 10~^° gcm"'^. Com- 
paring Figure [T] with Figure 2 of IMizuno et al.l (|l978l ) , the 
thermal evolution of the model with pcri = 10 pc,o (Fig. [T] 
broken line) corresponds to t hat for a gas envelo pe around a 
proto- Jovian planet (Fig.2 of lMizuno et al.l l l978f) when a gas 



opacit y Kg = 1.0 x 10 ^ cm^ g ^ is adopted. Mizuno et al.l 
l|l978t ) adopted Kg = 1.0 x 10"'' cm^ g ^ as the most reli- 
able parameter of a proto- Jovian planet, indicating that a 
more realistic gas temperature of the envelope is lower than 
that in the model with pcri = 10 pco (the dotted line of 
Fig. [1]). Note that the critical density pcri increases as the 
gas opacity Kg decreases. Thus, in models with pcri = lOpc.o, 



the thermal energy around the protoplanet may be overesti- 
mated. On the other hand, when the isothermal equation of 
state is adopted, the thermal energy around the protoplanet 
is obviously underestimated. Therefore, it is expected that 
the actual thermal evolution exists between models with 
Pcri = lOpco and pc = 00. 



2.3 Circumstellar Disk Model 



The initial settings arc similar to lMivoshi et al.l l|l999l) and 
[Machida et al. . (2006b , 2008, ). The gas flow has a constant 
shear in the i-direction as 



(0, --f^px, 0). 



For hydrostatic equilibrium, the density is given by 

2 



z 
'2h? 



(9) 



(10) 



where ao (= po dz) is the surface density of the unper- 
turbed disk. The scale height h is related to the sound speed 

Cs as h = Cs/f2p. 

In the standard solar nebular model (jHavashil Il98ll : 
iHavashi et al. |[l985l ). the temperature T, sound speed Cg, 
and gas density pc,o are given by 



-1/2 



(11) 



where L and are the protostellar and solar luminosities. 



kT 
prriH 



1/2 



= 1.9x10' 



/ T V''^ ( 2.34 

ViokJ \~ 



1/2 



where ^ = 2.34 is the mean molecular weight of the gas 
composed mainly of H2 and He, and 



Pco = 1.4 X 10 



-11/4 



gem 



(13) 



When Mc — 1 Mq and L = 1 Lq axe adopted, using equa- 
tions (|3]), (|lip . and (|12l) . the scale height h can be described 

as 

5/4 



/i = 5.0 X 10 



11 / Qp Y 
VlAU^ 



cm. 



The inverse of the angular velocity fip is described as 



0.16 



2/3 



(14) 



(15) 



The non-dimensional quantities are converted into dimen- 
sional quantities using equations (|liP " (|15p . 



2.4 Scaling 

The basic equations can be normalized by unit time, Slp^, 
and unit length, h. The density is also scalable in equa- 
tions (1) and (2) and is normalized by aa/h. Hereafter, the 
normalized quantities are expressed with a tilde on top, e.g., 
X — x/h, p = po / {(Jo / h) , t = tQp, etc. The non-dimensional 
unperturbed velocity and density are given by 

« = (0,-^5,0), (16) 



© 2008 RAS, MNRAS 000, [l]-?? 



4 M. N. Machida 



po 



2-K 



exp 



(17) 



Thus, non-dimensional equations corresponding to equa- 
tions JH), (O, and (0 are 



dp 
di 



-f V- (p*) =0, 



^ + (v ■ V)v = -ivP - V$cff -2z X V, 
dt p 

z r 



P^p 



1 _ tanh + l P 



tanh 



(18) 
(19) 
(20) 

, (21) 



where z is a unit vector directed to the z-axis. The gas flow is 
characterized by two parameters, the non-dimensional Hill 
radius fn = ru/h, and critical density pcri- In this study, 
the Hill radii fn = 0.29 — 1.36 are adopted. As a function 
of the orbital radius and the mass of the central star, the 
parameter ru is related to the actual mass of the protoplanet 
in units of Jovian mass Mj as 



Mo 



Mj 



: 0.12 



\Mr, 



Vl Al 



Auy 



(22) 



For example, in the model with fn = 1.0, ftp = 5.2 AU and 
Mc = 1 Mq, the protoplanet mass is Mp = 0.4Mj (model 
M04 in Table[T|. In the parameter range of rn = 0.29 — 1.36, 
at Jovian orbit (op = 5.2 AU), protoplanets have masses of 
0.01 — lAfj. Model names, non-dimensional Hill radii fn, 
masses of protoplanets at Jovian (5.2 AU) and Saturnian 
orbits (9.6 AU) and critical densities are listed in Table 1. 
Model names consist of two parts: the mass of the proto- 
planet at Jovian orbit and critical density. For example, 
model M001A3 has parameters of (Mp, pcri) = (O.OlMj, 
10^ Pco). 



3 NUMERICAL METHOD 
3.1 Nested-Grid Method 

To investigate the formation of a circumplanetary disk in a 
circumstellar or protoplanetary disk, it is necessary to cover 
a large dynamic range of sp atial scale. Using the nested- 
grid method (for details, see iMachida et al. I^OOsl . l2006al '). 
the regions near (using grids with higher spatial resolution) 
and remote from (using grids with lower spatial resolution) 
the protoplanet are covered. Each level of a rectangular grid 
has the same number of cells (= 64 x 128 x 16), but cell 
width As(/) depends on the grid level I. The cell width is 
reduced by 1/2 with increasing grid level [l I + 1). We 
use 8 grid levels (Zmax ~ 8). The box size of the coarsest 
grid, / = 1, is {Lx, Ly, Lz) = (12,24,3), and that of the 
finest grid, Z = 8, is (L^, Ly, L^) = (0.09375, 0.1875, 0.0234). 
The cell width in the coarsest grid, Z = 1, is As = 0.1875, 
and it decreases with As = 0.1875/2*"^ as the grid level I 
increases. Thus, the finest grid has As(8) ~ 1.46 x 10"^. The 
fixed boundary condition in the x- and z-direction, and the 
periodic boundary condition in the y-direction are adopted. 

The local simulation is not appropriate t o treat the gap 
formation because of the radial boundary l|Mivoshi et al. I 



ll999l : lTanigawa fc Watanabell2002l ). The gap property (the 
gap depth and width) depends on the size of the simulation 
box. Although the density gap may affect the mass accretion 
rate onto the protoplanet, it does not affect the formation of 
the circumplanetary disk, because the formation of circum- 
planetary disk depen d only on the size of Hill sphere (see, 
llMachida et al. II2008I ). 

3.2 Sink Cell and Smoothing Length 



In the finest grid (In 



the cell width is As — 



1.46 X 10~ . In real units, when the protoplanet is located 
at 5.2 AU, the cell width corresponds to As = 5.7 x 10^ cm, 
or 0.8 times the Jovian radius. In each model, the evolution 
of the protoplanetary system was calculated adopting the 
sink. The radius of the sink is fsink ~ 3.53 x 10~^ or twice 
the Jovian radius at Jovian orbit. During the calculation, 
the gas from t he region inside t he sin k radius is removed in 
each time step. lMachida et al. I (|2008l ) investigated the effect 
of the sink by parameterising the sink radius, and found 
that the circumplanetary disk and angular momentum of 
the protoplanetary system can be appropriately estimated 
when the sink radius is sufficiently smaller than the Hill ra- 
dius (r <C 0.1 rn). 

The smoothing length for the gravitational potential 
of the protoplanet is not explicitly adopted. In numer- 
ical settings, the physical quantities are defined at the 
cell centre, while the origin (protoplanet's position) is de- 
fined at the cell boundary. Thus, since the region inside 
r < fs = \/3As(Zinax)/2 has a uniform gravitational po- 
tential. At Jovian orbit, has a 0.7 times the Jovian radius 
(rs = 1.26 X 10"^ or = 4.9 x lO'' cm). 



4 RESULTS 

4.1 Models with Jovian-Mass Protoplanet 

In this subsection, to investigate the relationship between 
thermal effects and the structure and angular momentum 
of the protoplanetary system in detail, only models with 
a Jovian-mass protoplanet at Jovian orbit (models MlAl, 
M1A2, M1A3, Mil) are shown. 

4.1.1 Large Scale Structure and Density Distribution 

In each model, the evolution of the protoplanetary sys- 
tem was calculated for t ^ 600 (or JSlOO orbits), in which 
the flow around the protoplanet reached a steady state 
in a short timescale of t ~ 6 — 10 (~ 1 — 10 orbit). 
Previous st udies also show a steady state in a short 
times cale (e.g., iMivoshi et al. 1 ll999l : lTanigawa fc Watanabd 



l2002l : iMachida et al. Il2006d) . In this study, for safety, the 
evolution of the protoplanetary system was calculated for 
a sufficiently long time {t ^ 600) in all models. Since the 
evolutio n of the system until t he steady state was already 
shown in IMachida et al~l (|2008t ) in detail, only the evolution 
after the steady state is shown below. 

The left panel of Figure [2] shows the density (colour) 

^ The non-dimensional time unit i can be converted into the 
orbital period as f/(27r). 
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and velocity distribution (arrows) in the equatorial plane 
{z = 0) at t — 637.23 for model Mil, in which the gas be- 
haves isothermally. Model Mil has a mass of IMj at the 
Jovian orbit Op = 5.2 AU. To stress the structure inside the 
Hill radius, represented by the dashed line, only the grids 
of / = 3 — 6 are plotted in this panel. Note that the grid 
size of / = 3 is {Lx, Ly, L^) = (3, 6, 0.75), while the bound- 
ary of the outermost grid {I = 1) is located at {x, y, z) = 
(12,24,3). This panel shows the shocks (crowded contours 
near the Hill radius) in the upper right and lower left re- 
gion against the protoplanet (or the centre). Inside the Hill 
radius, a non-axisymmetric pattern appears. Those features 
are a l so seen in previous studies (K lcy 1999; Mivosli i et al. . 
,1999; Lu bow et anil999l : iKlev et al. 1120011 : iD'Angelo et al. I 
2002, . ,200^ ). 

The right panel of Figure [2] shows the density (colour) 
and velocity distribution (arrows) in the equatorial plane 
(z = 0) at i = 637.73 for model MlAl, in which the tem- 
perature increases gradually as the gas density increases, as 
shown in Figure [1] Figure [2] shows that the structures in 
models Mil and MlAl are nearly identical, while the shock 
in model MlAl appears in a more remote region from the 
protoplanet than in model MIL In addition, model MlAl 
has a more axisymmetric structure near the protoplanet 
(f < 0.5) than model MIL These differences are caused 
by thermal effects around the protoplanet. The larger ther- 
mal pressure forms a more spherical structure because the 
thermal pressure gradient force is isotropic. 

Figure 13] shows the density (colour) and velocity distri- 
bution (arrows) in the y = plane for models Mil, M1A3, 
M1A2 and MlAl, in which three levels of the grid (/ = 6, 7 
and 8) are superimposed. This figure indicates that models 
with higher gas temperature form relatively thick envelopes. 
The envelopes of models Mil and Ml A3 have a height-to- 
radius ratio of H/R ~ 1/10 in the region of p > 10^, while 
that of model MlAl has H/R ~ 1/5 in the region of p > 10^ 
(see, the contours for p = 10^ in Fig. However, even in 
model MlAl, the envelope in the region of p > 10" is suf- 
ficiently thin {H/R ~ 1/10). Thus, a thin disk forms close 
to the protoplanet in each model, while the protoplanetary 
systems have different envelope thicknesses in the region far 
from the protoplanet for models with different equations of 
state. In adiabatic models, since the gas temperature in- 
creases with the gas density as shown in Figure [T] the ther- 
mal effect is more significant close to the protoplanet. De- 
spite this, the structural difference is clearer in the region far 
from the protoplanet (i.e., the low-density region). This in- 
dicates that the protoplanet 's gravity and centrifugal force 
significantly dominate the thermal pressure gradient force 
close to the protoplanet. 



4.1.2 Angular momentum of Protoplanetary System 

Figure 13] shows the distribution of the average specific angu- 
lar momentum jV against the distance from the protoplanet 
f for models Mil, MlAl, M1A2 and M1A3. The average 
specific angular momentum jV is defined as 



Mr = / inr^pdr, 



and angular momentum 



4:TTr'^ pzui)^ df, 



(24) 



(25) 



Mr' 

where the mass 



(23) 



are integrated from the centre (f = 0) to a distance r. The 
solid line in Figure J?] shows the distribution of jV for the 
isotherm al model (Mill). This d istribution corresponds well 
to that of lMachida et alTl (|2008h . in which they investigated 
the accretion of angular momentum into a protoplanetary 
system using an isothermal equation of state without a sink 
cell, and found from the Jacobi energy, Kepler velocity and 
distribution of mass and angular momentum around the pro- 
toplanet that the angular momenta bound by protoplane- 
tary system are limited i n the region of f ^ 0.5fH — Irn. As 
shown in Figurellland Ma chida et al. I (|2008|). in each model, 
the average specific angular momentum increases from the 
centre to a peak around the Hill radius, and drops sharply 
just outside the Hill radius. The drop indicates that the an- 
gular momentum becomes negative at f > fn. Thus, the ro- 
tational direction is opp osite in region inside and outside the 
Hill ra dius. As shown in Sekiva et al. I (|l987f ). lMivoshi et al. I 
(|1999[ ). and iTanigawa fc Watanabd (|2002l i. a protoplanet 
formed by gas accretion in the circumstellar disk has a pro- 
grade spin, and thus it has a positive (specific) angular mo- 
mentum. On the other hand, gas far outside the Hill sphere 
seems to rotate retrogradely against the protoplanet because 
it rotates with nearly Keplerian velocity with respect to the 
central star [v = —(3/2) Qp x; see eq. As a result, gas in- 
side the Hill radius has a positive angular momentum, while 
that outside has a negative angular momentum against the 
protoplanet. 

Figure 13] shows little difference in average specific angu- 
lar momentum between the isothermal (Mil) and adiabatic 
(MlAl, M1A2, and M1A3) models in the region of f < 0.1, 
while there are no differences in the region of f > 0.1. Thus, 
the gas envelope near the protoplanet is slightly infiuenced 
by thermal effects, while differences in thermal energy do not 
affect the evolution of angular momentum in the region of 
f > 0.1. Since the angular momentum bound by the proto- 
planet is mainly distrib uted in the region of 0.5 fn ^ r Ifu 
(|Machida et al. II2008I ). the different equations of state (or 
different gas temperatures) hardly affect the total angular 
momentum acquir ed by a protoplanetar y system from the 
circumstellar disk. iMachida et al. I (|2008l ) noted that gravi- 
tational energy dominates thermal energy inside the Hill ra- 
dius when the mass of the protoplanet exceeds Mp > O.OSMj 
at Jovian orbit. Since models Mil, MlAl, M1A2 and M1A3 
have a mass of IMj at Jovian orbit, the gravitational energy 
dominates the thermal energy greatly. Thus, it is natural 
that different thermal energies hardly affect the evolution 
of the protoplanetary system and the acquisition process of 
angular momentum in these models. 



4.2 Relation Between Angular Momentum and 
the Protoplanet's Mass 



In i]4.1l we examined only the evolution of models having 
protoplanets with IMj. In this subsection, the evolution of 
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angular momentum for models with different protoplanetary 
masses is described. 

The average specific angular momenta derived from all 
models are plotted against the protoplanet's mass in Fig- 
ure (5] The average specific angular momenta are est i mated 
in the region of f < fii/2. Note that iMachida et al. I (|2008l ') 
showed that the angular momentum bound by the proto- 
planetary system is distributed in the range of l/2fH ^ 
f ^ f H and has almo st the same value in this range (for 
details, see Fig. 11 of iMachida et al. I [iool) . Figure [5] in- 
dicates that the system has almost the same average spe- 
cific angular momentum in models with the same proto- 
planetary mass but different equations of state when the 
mass of the protoplanet exceeds O.lMj. Thus, although 
the protoplanetary systems have different thermal energies 
around the protoplanet, they acquire almost the same spe- 
cific angular momentum from the circumstellar disk when 
Mp > O.lMj. On the other hand, when the protoplane- 
tary mass is smaller than O.lMj, the specific angular mo- 
menta for adiabatic models are slightly smaller than that for 
the isothermal model. Thus, the thermal pressure impedes 
the acquisition of angular momentum from the circumstellar 
disk when M < O.lMj. However, since a large fraction of gas 
and angular momentum falls into the protoplanetary system 
after the mass of the protoplanet becomes comparable to the 
present masses of gas giant planets, a difference in specific 
angular momentum in the lower-mass phase (A/ < O.lMj) 
hardly affects the subsequent evolution. 

The thin s olid l ine in Figure [S] denotes j oc M. 
IMachida et al. I (|2008l ) showed that, using the isothermal 
equation of state, the specific angular momentum is fitted 
by J oc M when the protoplanet has a mass of O.lA/j < 
M < IMj. Figure [5] shows that the average specific angular 
momentum can be fitted by j oc Af even wh en the adiabatic 
equat ion of state is adopted. As shown in IMachida et al. I 
(|2008l ). in real units, the average specific angular momen- 
tum of a protoplanetary system is approximated as 

when the protoplanet has a mass of Mp < lA/j. Although 
IMachida et al. I l|2008l ) estimated the angular momentum in 
the isothermal regime, equation (|26|l is valid even in the 
adi abatic reg i me as shown in Figure [5] 

iLissauerl l|l995h analytically estimated the specific an- 
gular momenturn 7'ana of a proto-planetary systern (see , also 
IStevenson et al. lllQSd and iMosgueira fc Estradall2003h as 

Jana = ^^^p?"!- (27) 

Using equations (|3]) and ([5]), equation (|27|l can be rewritten 
as 

'■'-'^"<">"{fkf imj)'" ^ ,2s, 

The different power of Mp and Op between equations (f26)) 
and (|28p is due to the disk model adopted. The standard 
disk model [see eqs. (|lip - (ll4|l ] is adopted in this study, 
while equation (|28|) is derived assuming only a Keplerian 
rotating disk in which the thermal pressure of the cir- 
cumstellar disk is ignored. Despite this difference, equa- 
tion (|28p is quantitatively the same as equation (I26|) . For 
example, with a proto- Jovian planet (Mp — lA/j) at Jo- 



vian orbit (ap=5.2AU), the protoplanetary system has a 
specific angular momentum of jim.j = 1-4 x lO^^cm^s"^ 
in equations p6|l . and jana,.T = 1.2 x lO^'^cm^s"^ in equa- 
tion (|27p . With a proto-Saturnian planet (A/p,j — 0.3A/j) 
at Saturnian orbit (ap=9.6AU), the specific angular mo- 
menta are jim.s ~ 1-2 x lO^'^cm^ s^^ in equation (|26p . and 
jana,.i = 7.2 X lO^'^cm^s"^ iu equation ^8^. Equations ((26)l - 
PSI) indicate that the specific angular momentum of a pro- 
toplanetary system is determined only by the protoplanet's 
mass (or the size of the Hill sphere) at a fixed orbit. Since 
the gravitational energy dominates the thermal energy in- 
side the Hill sphere when M > O. OSAfj at Jovian orbit as 
described in IMachida et al~l (|2008l ). the thermal effect can 
be ignored in the angular momentum acquisition process 
when the protoplanet's mass is Af > O.OSAfj. Thus, it is 
natural that the angular momentum fiowing onto the sys- 
tem is controlled only by the protoplanet's mass. 

The resulting angular momentum of the protoplanetary 
system can be estimated by integrating equation (|26p us- 
ing mass up to the present values for gas giant planets. As 
a result, the Jovian and Saturnian protoplanetary systems 
have ~ 30 and ~ 50 times larger angular momentum than 
the present gas giant planet syst ems (the central gas g iant 
planet + satellites, for details see IMachida et al~ll2008l ). 

4.3 Circumplanetary Disk around Proto-Jupiter 
and Saturn 

Figure [6] shows the density distribution on the y = plane 
(upper) and surface density along the z-axis (lower) around 
proto- Jupiter (left; model M1A2) and proto-Saturn (right; 
model M02A2). Model M1A2 (left panels) has a parameter 
of fu ~ 1.36, which corresponds to a Jovian mass (lAfj) 
at Jovian orbit (5.2 AU), while model M02A2 (right panels) 
has fu ~ 0.8, which corresponds to a Saturnian mass (IMs) 
at Saturnian orbit (9.6 AU). In addition, these models have 
Peri ~ lOOpcO, which is the most realistic param eter when a 
Jovian mass is adopted (see lMizuno et a"l]|l978l ). Since the 
evolution of protoplanetary systems was calculated in non- 
dimensional units, the calculation results can be scaled in 
real units at an arbitrary orbit. Figure |6] adopts the units at 
Jovian (left) and Saturnian (right) orbit. The dotted lines in 
Figure |6] represent 10, 30 and 50 times the planet's radius; 
a Jovian radius is adopted for model M1A2 (left panel), 
while a Saturnian radius is adopted for model M02A2 (right 
panel). 

The contours in the upper panels of Figure [6] show a 
thin disk close to the protoplanet (r ^ 50 rp), while the gas 
envelope far from the protoplanet (r > 50 rp) has a thick 
torus-like density distribution. In the region far from the 
protoplanet (r > 50 rp), the disk is flared, because the cir- 
cumplanetary disk connects smoothly to the circumstellar 
disk, which is thicker than the circumplanetary disk. Owing 
to rapid rotation near the protoplanet, the disk thins as it 
approaches the protoplanet. Figure [7] also shows the density 
distributions around the protoplanets in a bird's-eye view, 
but at a larger scale than Figure [S] Each surface in Fig- 
ure [7] represents an iso-density surface of p = 10'^ pc,o (red), 
100 PcO (orange) and lOpco (green and blue). The density 
distribution in the x — 0, y — and z = plane is projected 
onto each wall surface. Figure [7] shows that the central region 
sags in the centre of a concave structure, and a thin disk is 
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formed around the protoplanet. In addition, a butterfly- like 
structure is also seen on each wall surface. These structures 
are considered to be formed by the rapid rotation of the 
central circumplanetary disk. The thin disk near the pro- 
toplanet (red surface) is enclosed by a torus-like thick disk 
(orange and blue surfaces). Figures [S] and [7] indicate that 
a thin disk appears only in a compact region of r < 50 rp, 
and disk thickness becomes drastically large in the region of 
r > 50 rp. These features are considered to result from the 
angular momentum of the protoplanetary system. 

Using the centrifugal radius, the radial extent of the cir- 
cumplanetary (thin) disk can be related to the angular mo- 
mentum of the protoplanetary system. Under the assump- 
tion that the centrifugal force is balanced by the gravity of 
the protoplanet, the centrifugal radius is denoted as 



Jr 



(29) 



When the specific angular momentum and mass of the 
protoplanet Mp are given in equation (|29|l . the centrifu- 
gal radius can be estimated. Using the specific angular 
momentum of the protoplanetary system given by equa- 
tion (|26|) . the centrifugal radius for the proto- Jovian disk 
with Mp = IMj and ttp = 5.2 AU is 

rcf,j = 1.5 X 10" cm, (30) 

which is 21 times the Jovian radius (r = 21rjup, where 
rjup = 7.1 X 10® cm~^) or 0.028 times the Jovian Hill radius 
(r = 0.028 rn.jup, where rH,,iup = 5.3 x lO^^cm). In the 
same way, the centrifugal radius of the proto-Saturnian disk 
is estimated as 



rcf,s = 4.0 X 10 cm. 



(31) 



which is 66 times the Saturnian radius (r = 66rsat, rsat = 
6.0 X 10® cm) or 0.061 times the Saturnian Hill radius (r = 
0.061 rn.Sat, where rn.Sat = 6.5 x lO^^cm). These values are 
the same as the results of iMachida et alTI (|2008l ). in which 
the evolution of a protoplanetary system was investigated 
under the isothermal approximation with coarser spatial res- 
olution than this study. Observations showed that the regu- 
lar satellites around Jupiter and Saturn are distributed in a 
compact region in close proximity to the planet (r ^ 50 rp), 
and are on prograde orbits near the equatorial plane. The 
Galilean satellites, which are regular satellites of Jupiter, are 
distributed in the range of 6 rjup ^ r ^ 26 rjup, which is con- 
sistent with the centrifugal radius of equation (|30|l . The Sat- 
urnian representative regular satellites (Mimas, Enceladus, 
Tethys, Dione, Rhea, Titan, and lapetus) are distributed in 
the range of 3rsat ^ r ^ 60rsat, which corresponds well 
to equation (|3ip . In addition, the lower panels of Figure [6] 
show that the surface density has a peak in the region of 
lOrp < r < 50rp. This is as a result of gas accumulation in 
this annulus because of the centrifugal barrier. Thus, regular 
satellites may form in this compact region. 



5 CIRCUMPLANETARY DISK AND 
REGULAR SATELLITES 

After the formation of the circumplanetary gas disk, satel- 
lites may form in the disk. Altho ugh there are many sce - 
narios for satellite formation (e.g., IStevenson et al.lll986f ). 



regular satellites around gas giant planets are believed to 
form in the gas eous disk as planets form in the circum- 
stellar disk (e.g., ISteyenson et al. Ill986l : iKorvcanskv et alTI 
ll99ll : ICanup fc Wardll2002l ). Regular satellites are expected 
to form in the circumplanetary disk simultaneously with 
planets: the dust component condenses into a dust layer 
that fragments into satellitesimals in a thin disk near the 
proto planet, and these sat e llitesi mals merge to form satel- 
lites. iLunine fc StevensonI (|l982l ) suggested the minimum 
mass subnebular (MMSN) model, in which the mass of the 
present regular Jovian satellite system is reconstructed as- 
suming the dust-gas-ratio of the solar nebula. However, in 
this model, there are some difficulties for regular satellite 
formation. For example, planetesimals ~ 1 km in size fall 
into a central protoplanet wi thin ~ 1000 yr due to the gas 
drag (|Stevenson et al."l Il986l l. In addition, the disk itself 
is expect ed to be dissipated by viscous spreading within 
~ 10 vr ([Stevenson ct al. 1986). To overcome these difficul- 
ties, ICanup fc WardI l|2002l ) suggested the gas-starved disk 
model, in which the circumplanetary disk has several orders 
of magnitude lower surface density than in the MMSN disk 
model. In this scenario, solids delivered to the disk build 
up over many disk viscous cycles, resulting in a greatly re- 
duced gas-to-solids ratio during the final stages of satellite 
accretion. 

At the Callisto orbit (r ~ 26rjup), the MMSN dis k 
model has cr ~ 7.1 x 10 gcm~ (ILunine fc Stevenson|[l983 ). 
while the gas-starved disk model has a ~ 100 gcm"'^ 
(|Canup fc Wardll2002l '). Note that the gas-starved model as- 
sumes a very late stage of gas accretion onto Jupiter, in 
which the gas flows to the protoplanet very slowly. In this 
study, as shown in Figure[Sl the surface density at the Cellist 
orbit is (T ~ lO^gcm"^ which is an intermediate value be- 
tween those of the MMSN and gas-starved disk models. The 
Jovian representative regular satellites (Galilean satellites) 
are distributed in the region of 5.9 rp < r < 26.3 rp. This re- 
gion is sufficiently resolved in this calculation, since the cell 
width and sink radius are 0.8 rp and 2rp, respectively. How- 
ever, the inner boundary condition (i.e., the sink) may be 
inappropriate for investigating a circumplanetary disk close 
to the protoplanet. 

The sink cell is adopted to calculate the evolution of 
the protoplanetary system for a long duration. Without the 
sink, such calculation is impossible, because gas falls into the 
central region with a velocity much higher than the sound 
speed, and the time step becomes quite small. To properly 
estimate the circumplanetary disk in very close proximity to 
the protoplanet (r < 5rp), the protoplanet model (density 
and temperature) must be set at the centre. However, such 
treatment is beyond our scope: we cannot know the equation 
of state of the interior region of the protoplanet. Note that 
it is impossible to perform a calculation for the protoplanet 
itself, because it requires much greater spatial resolution. 

However, it is expected that the properties of the cir- 
cumplanetary disk do not change qualitati vely irrespective 
of the treatment of the inner boundary. iMachida et al. I 
(|2008t ) showed that the specific angular momentum and con- 
figuration of the disk are determined by the size of the Hill 
radius (or gravitational potential of the protoplanet). This 
study showed that the specific angular momentum of a pro- 
toplanetary system is independent of the thermal evolution 
around the protoplanet after the protoplanet's mass exceeds 
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O.lMj, indicating that the properties of the circumplane- 
tary disk are determined only by the outer boundary (i.e. 
the size of the Hill radius). On the other hand, the region in 
very close proximity to the protoplanet is slightly affected by 
thermal evolution, as shown in Figure U Since regular satel- 
lites are believed to form in this region, more advanced cal- 
culations are necessary to determine the circumstellar disk 
model closely related to satellite formation. 
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Table 1. Model parameters 



Model 




Mp* (5.2 AU) 


Mp* (9.6 AU) 


Peri 




MOOl (MOOlAl, M001A2, M001A3, MOOII) 


0.29 


0.01 


0.016 


(10, 100, 1000, 


oo 


M005 (M005A1, M005A2, M005A3, M005I) 


0.5 


0.05 


0.081 


(10, 100, 1000, 


oo 


MOl (MOlAl, M01A2, M01A3, MOII) 


0.63 


0.1 


0.16 


(10, 100, 1000, 


oo 


M02 {M02A1, M02A2, M02A3, M02I) 


0.8 


0.2 


0.35 


(10, 100, 1000, 


oo 


M04 {M04A1, M04A2, M04A3, M04I) 


1.0 


0.4 


0.65 


(10, 100, 1000, 


oo 


M06 {M06A1, M06A2, M06A3, M06I) 


1.15 


0.6 


1.0 


(10, 100, 1000, 


oo 


M08 {M08A1, M08A2, M08A3, M08I) 


1.26 


0.8 


1.3 


(10, 100, 1000, 


oo 


Ml (MlAl, M1A2, M1A3, Mil) 


1.36 


1 


1.7 


(10, 100, 1000, 


oo 



* in unit of Jupiter mass Afj 
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Figure 1. The density-temperature diagrams for isothermal 
(Pcri = oo) and adiabatic (pcri = 10pc,0i lOOpcOi and 1000pc,o) 
models. The right and upper axes indicate the dimensional tem- 
perature and density at the Jovian orbit, respectively. 
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Figure 2. Density {colour scale) and velocity distribution (ar- 
rows) on the cross section in the 2 = plane for models Mil (left) 
and MlAl (right). Level of the outermost grid is denoted in the 
upper left corner. Elapsed time i and velocity scale are given at 
the top of each panel. Dashed circle indicates the Hill radius. 
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\ = 637.23 .3 c logP s . ^ = 637.23 j , logp ^ 




Figure 3. Density {colour scale) and velocity distribution {ar- 
rows) on the cross section in the y = plane for models Mil (top 
left), M1A3 (top right), M1A2 (bottom left), and MlAl (bottom 
right). Level of the outermost grid is denoted in the upper left 
corner. Elapsed time i is given at the top of each panel. The ve- 
locity scale in units of the sound speed is denoted in the bottom 
of each panel. 
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Figure 4. Average specific angular momentum j against the 
distance f from the protoplanot for models Mil, MlAl, M1A2 
and M1A3. The vertical dotted line represents the Hill radius, 
f = 1.36. 
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Figure 5. Specific angular momenta for isothermal (+) and 
adiabatic models [p^ri = lOpc.o 100pc,o 10^pc,o (x) ] 

in the region of f < 0.5 rjj against r^, which corresponds to 
the protoplanet's mass [see eq. I|22|l ]. The upper and right axes 
indicate the dimensional mass in units of the Jovian mass and 
specific angular momentum at the Jovian orbit, respectively. 
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Circum-Jovian Disk Circum-Saturnian Disk 
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Figure 6. Circum-Jovian (left: model M1A2) and circum- 
Saturnian (right: model M02A2) disks around the protoplanet. 
The density distribution on the y = plane (upper) and surface 
density along the 2-axis (lower) are plotted. Dimensional units 
are at Jovian (5.2AU; left) and Saturnian (9.6 AU; right) orbits. 
The dotted lines represent 10, 30 and 50 times the protoplanet's 
radius, in which present Jovian (left panels) and Saturnian (right 
panels) radii are adopted. 
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Circum- Jovian Disk Circjm-Saturnian Disk 




2.0X10^ km 3.1x10^ km 

Figure 7. Structures around proto-Jupiter (left: model M1A2) 
and proto-Saturn (right: model M02A2) in a bird's-eye view. Den- 
sity distributions on x = 0, j/ = 0, and z = planes are plotted 
on each wall surface. The colour surfaces indicate iso-density sur- 
faces: p = lO^Pc.o (red), p = lO^PcO (orange), and p = lOpc.o 
(green and blue). Velocity vectors (arrows) are plotted on the 
bottom wall. Dimensional units are at Jovian (5.2AU; left) and 
Saturnian (9.6 AU; right) orbits. The size of the domain is shown 
in each panel. 
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